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' The motion of a particle in the field of dispiration (due to a wedge disclination and a screw dislocation) 

O |. is studied by path integration. By formulating a 50(2) x r(l) gauge theory, first, we derive the metric, 

curvature, and torsion of the medium of dispiration. Then we carry out explicitly path integration for 
the propagator of a particle moving in the non-Euclidean medium under the influence of a scalar potential 
and a vector potential. We obtain also the winding number representation of the propagator by taking 
the non-trivial topological structure of the medium into account. We extract the energy spectrum and the 
eigenfunctions from the propagator. Finally we make some remarks for special cases. Particularly, paying 
, attention to the difference between the result of the path integration and the solution of Schrodinger's equa- 

^ ' tion in the case of disclination, we suggest that the Schrodinger equation may have to be modified by a 

\l , curvature term. 

, PACS: 61.72Lk - Quantum mechanics, linear defects: dislocations, disclinations. 
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1 Introduction 



In recent years quantum effects on particle propagation in a field of topological defects have attracted 
I considerable attention (see, e.g., [U El E]; see also a recent review article [4]). Although the notion of 

"defects" in physics was originally associated with crystalline irregularities, it has been extended to more 
general topological structures such as entangled polymers, liquid crystals, vortices, anyons, global monopoles, 
cosmic strings, domain walls, etc. In 1950's, from the structuralogical aspect, Kondo et aZ. [5] extensively 
studied unified geometrical treatments of various subjects including elastic and plastic media, relativity, 
network systems, and information theory. In particular, Kondo [6] related dislocations to Cartan's torsion 
in the medium. Since then the relation between dislocation theory and non-Riemannian geometry has been 
well established. More recent approaches to defect problems are gauge theories similar to those in particle 
physics [71 [H m [To]. In 1978, Kawamura [11] pointed out that a screw dislocation in a crystal produces 
an Aharonov-Bohm type effect in particle scattering. The Aharonov-Bohm effect [12] is usually understood 
as a topological effect [TS] [HI [151 US] • Since standard approaches to particle-defect interaction problems 
constitute of solving relevant (local) differential equations, the role of topology becomes often obscure. 
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In the present paper, we analyze the quantum behavior of a particle in the vicinity of a topological 
defect by thepath integral method. Specifically we carry out path integration for a particle moving around 
a dispiratiorp. The dispiration we consider is a composite structure of a screw dislocation and a wedge 
disclination with a common defect line. In Sec. 2, we study the geometrical and topological properties of 
the dispiration by using the 5*0(2) x r(l) gauge theory. For the medium M of the dispiration, we derive 
the squared line element, 

ds^ ^dr^ +a^r^d0^ + {dz + pd0f, (1) 

where /3 and a are the parameters directly related to the Burgers vector of dislocation and the Frank vector 
of disclination, respectively. The medium has a non-Euclidean structure with singular torsion and curvature 
along the dispiration line. 

In Sec. 3, we calculate in detail the path integral to obtain the propagator (Feynman kernel) for a particle 
in the field of the dispiration characterized by the line element ([1]) under the influence of a scalar potential 
and a vector potential. From the propagator so obtained, we extract the energy spectrum and the energy 
eigenfunctions. In Sec. 4, converting the propagator in the partial wave expansion into that in the winding 
number representation, we show that the effect due to the multiply connected structure of the medium is 
taken into account. The energy eigenfunctions and the energy spectrum are extracted from the propagator 
in Sec. 5. Sec. 6 is devoted to interpretation of the energy spectrum. In particular, we discuss the difference 
between the results from the Schrodinger equation with no curvature term and those from the path integral 
for the case of the conical space. Our path integral calculation suggests that the standard Schrodinger 
equation may have to be modified by a curvature term in order for two approaches to be consistent. 



2 Gauge formulation of dispiration 

The dispiration under consideration is a combined structure of a screw dislocation along the z-axis and a 
wedge disclination about the same z-axis. See Figs. 1-2. The gauge theoretical approach to dislocation 
and disclination has been extensively discussed in the literature [II [H [9l [10] . Here we wish to present a 
SO{2) X r(l) gauge approach to the dispiration developed along the line similar to the formulation by 
Pmitigam and Soleng [lOj . 

The S0{2) X r(l) gauge transformation:— In the gauge theoretical treatment, a deformation of an elastic 
medium in three dimensions is described by a local coordinate transformation consisting of a rotation and a 
translation, 

x' = p(x)-x + r(x), (2) 

where p G S0{2), t G 2^(1), and x ~ {x,?/, z} is the position vector of the undeformcd three-dimensional 
medium. 

As is well-known, an axial wedge disclination can be created by the so-called Volterra process: That is, 
by (a) removing a wedge shaped portion of an angle 7 G [0, 27r) from the medium and pasting the open walls 
together, or (b) inserting an extra wedge-shaped portion with 7 G [— 27r,0) into the medium. Note that the 
deficit angle 7 is chosen to be positive for the case (a) and negative for the case (b) . See Fig. 1 for the case 
(a). The wedge disclination about the z-axis is a rotational deformation obtained by gauging 50(2): 

x' = p{Q) ■ X (3) 

with the rotation matrix, 

cos(76'/27r) - sin(76'/27r) 
p{e) = \ sin(76i/27r) cos(76l/27r) | . (4) 
1 



^The Volterra process of forming a dislocation and a disclination involves a translation and a rotation, respectively. 
A dispiration is a defect formed by a translation and a rotation at the same time [17]. See Fig. 1 and Fig. 2. 
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Here = ta.n-^{y/x) G [0, 2tt). 

A screw dislocation lying along the z-axis (i.e., having a constant Burgers vector pointing the z-direction) 
as shown in Fig. 2 is a translational deformation obtained by gauging the z-translational group T{1): 

x' = x + t(0) (5) 

with the angle-dependent z-translation vector 

rid) = -^e, (6) 

where 6 is a translation parameter and 9 = tan^^(?//a;) as before. 

The dispiration comprised of such a wedge disclination and a screw dislocation is described by the 
combined coordinate transformation, 

= p{e)-x + T{e) (7) 

where p{9) and t(6') are given by ^ and ([6]), respectively. 

Gauge connections:— The standard Yang-Mills theory localizes with respect to the external space the 
gauge group which acts homogeneously in the internal space and introduces the gauge potential or the 
connection T that transforms under the group action Q as 

T' = grg-^ - dg g-\ 

What we wish to formulate for the dispiration is a 50(2) x r(l) gauge approach which differs in character 
from the Yang-Mills theory; the gauge group acts inhomogeneously on the (internal) coordinates x that is 
soldered locally to the (external) coordinates of the medium of dispiration. It is more appropriate for us to 
follow the procedures used in constructing a Poincare gauge theory [lUl [TB] or an afhne gauge theory [T^] for 
gravity. 

First we write the group action ([7|) in the matrix form. 



^x- 

Then we define the connection T as 







where F^^^ and F^^^ are the rotational connection and the translational connection, respectively. In order 
for F to behave like a connection in the Yang-Mills theory, F^^^ and F^-^^ must transform as 

F(^)' = pF(«)p-i-(dp)p-i 

and 



,-1 



F(^)' = pF(^) -dT- pF(«)p-i - (dp) p 

Evidently the rotational connection is a 50(2)-valued one-form, and the translational connection is a R- 
valued connection one-form. 

Now we construct the solder form that locally connects the global gauge coordinates x(0) to the coordi- 
nates x(0) of the dispiration medium as 

= dx + F(^) • X + F(^) (8) 
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which is a vector valued one-forni0 It transforms as oj' = pu3, leaving g — uP^ ■ u} invariant. Since we create 
the dispiration by a gauge transformation in flat space, we choose a connection that vanishes at — 0. Then 
the corresponding rotation and translation connections are given, respectively, by 

r(^) = pdp-i (9) 

and 

r(^) = ~dT. (10) 

The rotational connection ([9]) and the translational connection ([T0| can be easily calculated by using the 
rotation matrix (j4|) and the translation vector (|6]). The differential of the rotation matrix (jlj is 



27r 

where 



sin(76'/27r) - cos(76'/27r) 
76'/27r) sin(76'/ 




dp-'^{e)^dp{-e)^-^de \ cos(76i/27r) sin(76'/27r) I (11) 



de =^ \{xdy - ydx) , = + (12) 



This and the rotation matrix ([4]) together lead us to the rotational connection, 



7 

± ■ - = py^a) up yu) = 

where 



T'^^^ ^ p{0)dp-\0) ^ -^mdO (13) 
27r 



1 

m = I - 1 I . (14) 


The translational connection is found in a simple form by differentiating the translation vector 

T^^^ ^ -dr ^ ^d9e,. (15) 



Substitution of (|T3l) and (|T5j) into dS]) yields 
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, dx + {-f/2Tr)yd0 

dyi+ —ddm-x+ —d9e^^ \ dy - l-//2Tr)xd0 \. (16) 

\dz + {b/27r)de 



Utilizing this solder form as the coframe we can determine the squared line element ds^ = gij dx'' (8> dx-' in 
the medium of the dispiration, 

ds^ = ScfiUj" (S)UJ^ ^ dr^ + d0^ + {dz + I3d0f , (17) 

where 

Evidently a e (0, 1] when 7 e [0, 27r), and a e (1, 2] when 7 e [-27r, 0). 



■^It is possible to derive this solder form by starting with the homogeneous group 50(4) and by applying the group 
contraction to reduce it to SO{3) x r(l). The solder form may be obtained as a limiting form of a component of 
the connection. For instance, the de Sitter gauge theory can be contracted to the Poincare gauge theory, so that the 
solder form is obtained from the gauge potential as the vierbeine [20) . This contraction scheme does not work for the 
afiine gauge theory |19) . 
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In the present paper, we consider the non-simply connected medium — R"^\{a; = y — 0} with metric 
(fTT]) as the field of dispiration. In the following, we briefly review some of the geometrical and topological 
properties that will be useful for later discussion. 

Curvature and Frank vector:— The curvature two-form R is defined in terms of the rotational connection 

R = dr^-^) + r(^) Ar(^). (19) 

For the case of dispiration with ([13]), since r*-^-* A r*-^-* = 0, it integrates into 



[n^fdT^^^^f r^-^) = / -^md6l = 7m, (20) 
Js Js Jas Jas ^tj" 

where dS denotes the boundary of a surface S. Now we choose an orthogonal frame such that 

ds^ = Ski^'' A with ^1 = dr, ^2 ^ arde, ^ dz + pdO. Then we have R = Rfei^'' A The surface is 
chosen to be orthogonal to Eq. (PH)) implies that 

R = 7m5(2)(^i^^2^rf^i ^21) 

The corresponding scalar curvature is 

R = 2^5^^\e.e)- (22) 

or equivalently 

R = 2l5(^\x,y). (23) 
a 

Apparently the curvature of the medium is zero everywhere except along the z-axis {x = y = Q). Note 
that the curvature is created not by the dislocation but by the disclination. 

The Frank vector f = {$^'^, $^^} is defined to characterize a disclination with 



y R = 7 m (24) 



where again the integral is over a surface S which is delimited by a loop dS enclosing the z-axis. Obviously 
the only non- vanishing components of the Frank vector for the dispiration is given by As a result, the 
Frank vector takes the form, 

f-7e„ (25) 

which, points the z-direction and its magnitude is identical to the deficit angle 7. 
Torsion and Burgers vector:— The torsion two-form is defined by 

T^du} + T^"'^ A w = R • X + dr^^' + r^-^' A r^"^^ , (26) 

which, in general, depends not only on the translational connection but also on the rotational connection. 
For ([13]) and ([T5]), T^^^ A T'^^ = 0. Note also that the integration of R • x vanishes. Then the torsion 
two-form integrates into 

' dr(^) = / r(^) = be,. (27) 
from which follows 

T = bej'^'^\£},i^)d£} hd^'^. (28) 

We see that the only non- vanishing component of the torsion two- form is in the z-direction and is contributed 
only by the screw dislocation. 



5 



The Burgers vector b is defined by the surface integral of the torsion two-form 

(29) 



which has been evaluated in (|27l) 

h^be^. (30) 
As is expected the translation parameter 6 of ([6]) is indeed the magnitude of the Burgers vector b. 

Conical space:— Before closing this section, we consider a special case of the line element ([1]) with constant 
z and /3 = 0. The two dimensional surface having the metric 

df = dr^ + a^r^de^ (31) 

to which ([T]) reduces may be realized as a conical surface Mc if the surface is imbedded into a three dimensional 
Euclidean space E^. Let 

X = (err cos d, ar sin 6*, v^l - 0-2 r) . (32) 



Apparently, 

I : dX-dX = dl'^ (33) 

which is the first fundamental form of the imbedded surface Ale- Again, parameterizing the surface by 
< r < cx) and < 9 < 2tt, we have the metric tensor and its inverse, 

9a^-(l I2 V ff"" = f J ] , (34) 



which will be used later. With 

Xr = {a cos 0, a sin 0, \/ \ — , (35) 

Xe = (—err sin 6',CTr cos 0,0) , (36) 
the unit vector normal to M^, at a point (r, 9) , 

n = Xr X Xe/|X^ X Xe| (37) 

is easily calculated to be 

n= (^-Vl -cr2 cose*,- Vl-cr2 sin 0,0-) (38) 

The second fundamental form, 

// : - dX • dn = Gabdu" ® d^t^ (39) 
is also immediately obtained for the conical surface Mc as 



- dX • dn = y/l - (t2 err ^6*2. (40) 

Hence 







Vl-cr2/(crr) 



(41) 



whose two eigenvalues fci = and k2 = vT—a^/ {ar) are the principal curvatures of the conical surface for 
r 0. The Gaussian curvature and the mean curvature of a two dimensional surface imbedded in are 
defined, respectively, by 

K = detG^kik2 (42) 
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H=\trG^\{k,+k2). (43) 



For the conical surface in question, 



K^O and H = ^Jl- a'^/{2ar) (44) 

for r 7^ 0. Of course the Gaussian curvature K does not vanish at the apex of the cone. According to the 
Gauss-Bonnet theorem, 

/ Kda + / kgdl = 2Trx{S) (45) 

Js JdS 

where S" is a compact two-dimensional Riemann manifold with boundary dS, K is the Gaussian curvature 
of S, kg is the geodesic curvature of dS, and x{S) is the Euler characteristic of S. Now we let S be the 
lateral surface of a frustum with slant height r. Then kg = 1/r and the line element dl integrated along the 
boundary dS with dr ~ result in 

/ kgdl = 2n(7 = 27r - 7 (46) 
Jas 

where 7 is the deficit angle defined in (jlSp . The Euler characteristic of a cone is unity. Hence we have 

Kda = 7 (47) 

JS 

where da = ardr dO. From this and (|44)) follows 

K=l5^(x,y) (48) 
a 

which includes the curvature at the apex {r — Q). As is well-known, the Ricci or scalar curvature of a two- 
dimensional manifold is twice the Gaussian curvature. Indeed, twice the Gaussian curvature (|48p coincides 
with the scalar curvature (1231). 



3 Path integration 

The medium surrounding the dispiration considered in the preceding section is = R'^\{x = y = 0}, which 
is geometrically a non-Riemannian space and topologically a non-simply connected space. In this section, 
we carry out path integration for the propagator of a particle moving in the field of the dispiration. To 
confine the particle in the vicinity of the dispiration, we introduce a two-dimensional harmonic oscillator 
potential. Furthermore, we assume a repulsive inverse-square potential to prevent the particle from falling 
into the singularity at the defect line. For the purpose of comparison, we introduce also a vector potential 
due to a flux tube and a uniform magnetic field. 

The standard approach deals with the Schrodinger equation in curved space. As will be discussed in Sec. 
5, the Schrodinger equation is usually modified in curved space not only by the Laplace-Beltrami operator 
replacing the Laplacian but also the so-called curvature term added as an effective potential. The energy 
spectrum is sensitive to the type of curvature term; yet the controversy on the choice of the term is not fully 
settled. The path integral calculation we present suggests that the Gaussian curvature of the surface where 
the particle moves would dominate the curvature term. 

3.1 The Lagrangian 

Now that the dispiration field is characterized by the line element (fT7|) . the Lagrangian for a charged point 
particle of mass M moving in the vicinity of the dispiration under the infiuence of a scalar potential ^(x) 
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and a vector potential A(x) is written as 

L = -M 



-i- A(x)-F(x). 



(49) 



As has been mentioned before, we choose the vector potential A(x) consisting of two parts; one due to an 
ideally thin flux tube that contains constant magnetic flux $ along the dispiration line, and another due to 
a uniform constant magnetic field B = _B 62 pointing the z-direction. Then the vector potential term of the 
Lagrangian (|49p is expressed in cylindrical coordinates (r, 0, z) as 



-X • A = aM + Muj^r'^L 
c 



where 



6$ 



eB 



(50) 

2Tihc 2Mc ^^'^^ 

Here a is the ratio of the magnetic flux to the fundamental fluxon $0 = 27r?ic/e, which is identical to the 
statistical parameter in Wilczek's anyon model |35j , and is the Lamor frequency. The flux tube is included 
in our calculation in order to observe the similarity between the Aharonov-Bohm effect and the effect of a 
screw dislocation. In the scalar potential T^(x), we include a two-dimensional short range repulsive potential 
(the inverse square potential with k > sufficiently large) to emphasize the impenetrable feature of the 
central singularity, and a two-dimensional long range attractive potential (the harmonic oscillator potential) 
to confine the particle in the vicinity of the dispiration. Namely, V{n) is specified to be a two-dimensional 
central force potential, 



2 ° 



y 



Thus the Lagrangian we consider is 

I 
2 



- afiB - Muj^ r^e ~ V{r). 



(52) 



(53) 



3.2 The propagator 



The transition amplitude (propagator) for the three-dimensional motion of the charged particle from point 



x' = (r', 6', z') to point x" = (7 



,// all ^11 



z") can be calculated by the path integral [2T] 



X(x",x';t) 



|.x"=x(t") 






/ exp 


- /* Ldt 




/x'=x(t') 


Tt- Jt' 





(54) 



where t = t" — t' > 0. The integral measure must be so chosen that the propagator satisfies the properties, 

lim K (x", x'; t" - t') = 5(x" - x'). 



/A-,x".x;<"-„A-,x.x'.<-<V=x^i.-(x",x';, 

The path integral we calculate with the Lagrangian ([55)) is 
K{r",z";r',z';T) 



0. 



(55) 
(56) 



exp 



t" 



ahe - Muj^e ~V{r) \dt 



V^yVz. (57) 



In (|57p r signifies two variables (r, 9) symbolically, and the two dimensional integral measure (Pr will be 
specified later. After path integration, from the propagator, we should be able to extract the energy spectrum 
and the wave functions for the system. 
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3.3 The ^-integration 

First we perform the z-integration by letting — z + /39. The z-path integral is nothing but the Gaussian 
path integral for a one dimensional free particle, which yields the standard result, 



K"=C(t") 




/ exp 









M 
2TrihT 



exp 



iM ic" - C'f 

2n^ 



Now we rewrite the right hand side of ([55t as 

'iM{C' -Cf 



M 

2TTihT 



exp 



2hT 



1 

27r 



where hk is the z-componcnt of momentum of the particle. We also notice that 



C" - C' = z" -z' + P' {9" - 9') = z"-z' + (3 



9dt. 



Incorporating these results into the path integral (j57l) , we decompose the propagator as 

X(x",x';t) = — / dfce*'=(""-"')e-'"'"''''/2*''x('=) (r",r';T) 

with the two-dimensional propagator for a fixed k value, 
ifC^) (r",r';T) 



(58) 



(59) 



(60) 



(61) 



-"=r(i") 



exp < - 



r'=r(t') 



t" 



n.L, 



— {r^ + a^r^9'^ ) - Muj^r^9 - S,h9 - V{r) 



dt \ V^Y. 



(62) 



where ^ = a — pk. The integral on the right hand side of (j59p is the spectral decomposition of the z-motion 
with the continuous spectrum, 

Ek = (63) 



2M 

Next we make a change of the angular variable from to by letting 

d = 9-uj 

where Cj — uj^ ja^ . Accordingly the /c-propagator (|62p is transformed into 
ifW (r",r';T) 



(64) 



•"=r(t") 



exp ■ 



r'=r(t') 



I 



t" 



M 



P + a^r^d-" ] - iM- U{r) 



dt ) V^r{r, d) 



U{r) = V{r) + -MCj^r^ - ^fiu), 



(65) 

(66) 
(67) 
(68) 

Note that "d is the angular variable in a rotating frame with angular velocity —uj. For simplicity the constant 
Vq in the potential will be ignored in the calculation below. 



where 



or 



where 



U{r) 



1 



SMo-V^ 2 



uj'^ = uj^ + cj^, Vq — —S^huj. 
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3.4 The short time action 

To calculate the two-dimensional path integral (|65l) in polar coordinates [22j [231 El] , we first express it in 
discretized form, 

/^W(r",r';r)= hm / 1[K^'^^ {r,,r,.,;e) H d'r, (69) 

where the propagator for a short time interval e — tj — tj-i — t/N is given by 

i^W(r,,r,_i;e)=A,exp(^j5,). (70) 

Now we select a relevant approximation of the short time action Sj in cylindrical coordinates, and determine 
the amplitude Aj so as to meet the normalization condition 

hm ifW (r„r,_i; 6) = <5(2)( _ ) (7^) 
where the two-dimensional delta function satisfies 

J(^)(r,-r,_i)d^r, = l. (72) 



The short time action is 

r A/f . .... . 1 

dt (73) 



which we approximate by 



Sj^^{ {Ar,f + 2a\,r,^^ [1 - cos(A^,)] } - ^HA^, - U,e , (74) 

where Avj = Vj — rj-i, At9j = -dj — dj-i, and Uj — U{rj). It is tempting to approximate [ddY by (Az9)^. 
In path integration, however, {A§Y /e cannot be ignored because (Ai?)^ ~ e. Hence it is appropriate to 
replace (di?)^ by 1 — cos(A'!9) even though unimportant higher order terms are included (see [25]). Since it 
is sufficient for a short time action to consider the contributions up to first order in e, we further employ the 
approximate relation for small e |15| , 

cos (Ai?) -t- aeAi? ~ cos (Ai9 - ae) + ^a^e^ (75) 

to write the short time action multiplied by {i/K) as 

We use this short time action (multiplied by i/fi) for evaluating the angular path integral. However, before 
starting the angular integration, let us determine the amplitude Aj by considering the limit e — >■ of the 
short time action (l76l) . 



limA.e'^^/'' = limA,exp|^(Ar,)2| exp. 



A,J^S{Ar,)5{A^,). (77) 
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Now let the areal element be given by 

(frj = Tj drj ddj. 
Then the amplitude meeting the condition (|7ip must be of the form, 

Ma 



A, 



2T:ihe 



Alternatively, if 

then the corresponding amplitude is 



(78) 

(79) 
(80) 

Since path integration with either combination yields the same result, we employ the first choice (j78p with 
(17^ for our calculation. 

3.5 Asymptotic recombination 

For path integration in polar coordinates, separation of variables is not straightforward. To separate the 
angular variable from the radial function, we employ the asymptotic recombination technique (see [25]). The 
asymptotic form of the modified Bessel function I,^{z) for large \z\ (GR:8.451.5 in j26j) is 



(frj = arj drj ddj, 



M 



'2ttz 



> (2z)" ^ 



where —'in/2 < arg z < 7r/2, and 



(i^, n) 



V{v + n+\) 



12-KZ 



(i.,0) = l. 



(2z) 



(82) 



The asymptotic recombination technique is based on the conjecture that the one-term asymptotic form, 

1 



12-KZ 



exp 



(83) 



is valid for sufficiently large \z\ and for — 7r/2 < arg z < 7r/2 as relevant in path integration [221 127) . With 
the help of the one-term form which we refer to as the Edwards- Gulyaev asymptotic formula, we can derive 
the following asymptotic relation for large \z\ [25] . 



6z -c/2 



lv{az) e e 

where a > 0, 5 > 0, 8a(a + 6)c > 5 - 4(a + h)v'^ and 



a + 6 



[(a + h)z\ 



a + b 



h 2(a + h)c 

a 4a 



1/2 



Making use of this asymptotic relation together with the Jacobi- Anger expansion formula 

oo 



(84) 



(85) 



(86) 



m— — oo 
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we obtain another asymptotic relation for large z, 



exp 6z + z cos 



d 



M + i- 



(87) 



m— — oo 



with 



/i= [(l + 6){(m + d)2 + 2/}-fe/4]'/' 
In the above we have let a = 1 and c — {(P +2f + 2md)/2. 



3.6 Angular integration 

Utilizing (l87t with b — a^^ — 1, z — M(T'^rjrj^i/{ihe), d — ^ — a ~ f3k and / = k/2, we separate variables 
of the short time propagator ([70l) as 

i^W (r„r,_i;e) = A, exp (^S^) = £ 



m, — — oo 



with the short time radial propagator, 

M 



Rm, {fj,rj-i\e) = — exp 
ihe 



iM 
We 



4n 



Himj ) 



ihe 



where 



1 



fiimj) = — [4(toj + + cr^ - 1 + k 



1/2 



(89) 



(90) 



(91) 



The prefactor \/l + b appeared on the right hand side of the formula (f87|) results in a multiple factor 
before the summation of ([M]). which will cancel out the factor cr appearing in the amplitude (|79l) . Thus the 
two-dimensional path integral (1691) is reduced to the form, 



The angular integration can be done straightforwardly by using the orthogonality relation. 



(92) 



(93) 



^ira{'d —'d ) 



(94) 



Namely, 

JV-l N-l 
]-J^.™,A^, -Q ^ (2^)Ar-l -Q ^^^^^ 

where m = ttin- After angular integration, we have rrij — m for all j, and arrive at the expression for the 
full propagator with a fixed wave number fc. 



oo 

K^^){v",v'-t) = — V e,,v[im{r -d')]R^{T"y-T) 

ZTT ^ — ' 



(95) 



„ N N-l 

Rm{r", r'; r) = lim / J]^ i?™ (r^, r^-i; e) J| r^drj 

'J A 1 A 1 



where 

AT AT-l 

(96) 

The left hand side of (jOSp is nothing but the partial wave expansion in two dimensions, and the radial 
propagator ([M)) corresponds to the m-th partial wave propagator. However, the last radial path integration 
(|96p remains to be carried out. 
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3.7 Radial path integration 



To perform the radial path integration exphcitly, we first notice that the short time radial propagator (|90 
with m = ruj for all j is similar in form as that of the radial harmonic oscillator. Then we rewrite (|96p as 



M 

R„i[rj,rj^i;e) = — exp 
ine 



^MUJ 2 , 2 N 1 1 2 2 



fj.{rn) 



\ ihe 



We further simplify (^7]) by letting r]j = {Muj/2h)r^ and (p = arcsin(a;e) as 



i?„i(rj,rj_i;e) 



MuJ 

ih sin if 

Here we have used the approximation, 



exp [i{r]j + Tij-i) cot ip] (-2i^/r]Jr]JZ^ csap) 



22 



cosifj = cos[arcsin(aje)] ^ 1 — —uj e 



(97) 



(98) 



(99) 



At this point, for convenience, we introduce a two point function, referred to as the u-function |25[ 128] . by 



{rj, ri']ip) = —i CSC ip exp [i (77 + 77') cot ip] (^—2i^/r]rf esc ip^ 
satisfying the convolution relation 

poo 

/ iv", m f) t^/x (??, f) d-q = Vf, (77", 77'; 2ip) 
Jo 

which can be derived from Weber's formula (GR: 6.633.2 in [26]) as modified in [23] . 







exp(zar^) /p(— iar) /^(— i6r) r — — exp 



I {- — 
I 2a 



(100) 



(101) 



(102) 



valid for Re a > and Re /i > —1. Then the short time radial function can be expressed in terms of the 
u-function as 

Muj 

Rmirj,rj-i;e) = — {r]j,r]j-i;<p) . 
Now we are ready to perform the radial path integration. Substitution of this into (j96p gives 

N N-l 



(103) 



Rra{r",r';T) 



Muj 



The convolution property (|10ip of the ?;-function enables us to reach 

Era {r'\ r'; r) = ^^t^,,(m) (?/',?/; wr) . 
In the above we have also used the property, 

lim (Np) = lim (N sin p) = lim (ujNe) = lot. 

N—^oo N^oo N—>-oo 



In terms of the modified Bessel function it is written as 



Rm {r",r'-T) 



Muj 
i?isina; 



■ exp 



iMuj 
2h 



(r'2 _|_ j."2-) cot UT 



Mujr'r" 
ih sin wi 



(104) 



(105) 



(106) 



(107) 
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with 

^{m) = -!-v/4(to + a- likf + 0-2 - 1 + k. (108) 
2ct 

Notice that /i(m) is a real positive number \i 1 — < k. In this raanner we have completed the radial path 
integration for the partial propagator with k fixed. It turns out that the radial propagator we have obtained 
above is identical in form with that of the radial harmonic oscillator. The characteristics of the dispiration, 
the flux tube, the uniform magnetic field and the assumed potential, which make the present system different 
from the simple harmonic oscillator, are all taken into the index /i(m) of the modified Bessel function. 
With the radial propagator (jl07l) , the full propagator with k fixed is obtained by 

i^W (r",7?";r',7?';T) 
Mijj 

= - — — ex 

ZTTiasmwr 

or by converting § into 9 = "d + ut, 

i^W ir",e";r',e';T) 

Mlu fiA 
= — ; exp — 

2'Kinsm.ujT 2 



iMuj 

2h 



^) C0t( 



m— — 00 



yihsiiiujT 



(109) 



(r'2 + r"2) cotwT 



00 

E 

m— — 00 



,im{e" -e' -ujt) t 

fi\m) 



Mujr'r" 
iTismujT 



(110) 



4 Winding number expansion 

In the angular path integration performed above, we have not explicitly taken account of the topological 
structure of the background medium M.. Since M = R'^\{x = y = 0} = R x R+ x 5*^, the paths connecting 
two points, say a and 6, in Ai can wind around the z-axis many times, and may be classified (into homotopy 
classes) by the fundamental group 7ri(A^) = 7ri(S'^) = 7i oi M. A set of all homotopically equivalent paths 
in M. is now characterized by a single winding number n G Z. Therefore the propagator K{b,a) may be 
given as a sum of subpropagators Kn{b, a) for the paths with different winding numbers n, 

Kib,a) = J2CnKn{b,a). (Ill) 

TlSZ 

The paths in a same class may be deformed into one another, so that the transition amplitudes corresponding 
to these paths must share the same phase factor, On the other hand, the paths belonging to different classes 
may have different phases. Laidlaw and deWitt [29], and Schulman \TT argue that the coefficients C„ are the 
one-dimensional unitary representations of the fundamental group 7ri(Al) = Z. From cq. (jllip it is apparent 
that Cn+m — CmCn- If uo degeneracics and no internal degrees of freedom are assumed, Cn+i — e^^Cn- 
Then the coefficients are given by the one-dimensional representations of the fundamental group T:i{Ai) or 
the additive group Z, namely C„ = e*"''. 

It has also been pointed out |30) that the angular momentum representation and the winding number 
representation are complementary to each other via Poisson's sum formula, 

riGZ m6Z 

This means that the propagator K'^^\y" , r'; r) we have obtained in the preceding sections, as is given in the 
angular momentum representation, cannot be regarded as a subpropagator carrying a single winding number. 
Similarly, any partial propagator with a fixed angular momentum quantum number cannot be decomposed 
to a full set of subpropagators with winding numbers. The propagator obtained for the dispiration field must 
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be understood as a full bound state propagator that contains contributions from all homotopically possible 
paths (see [T51I16] for detail). Thus we look for a winding number representation of K^''\r" ,r']T) in ()109|) 
via Poisson's formula ()112|) . 
Let us rewrite (^5]) as 

i^W(r",r';r) = ^ f + A - m) e^("'+^)(''"-^') i?„,+,(r", r'; r) dA (113) 

'' mez 

with the radial propagator Rm{r" ,r' ;t) given by (|107|) . Then we utilize Poisson's formula (|112l) to convert 
(|113p into the winding number representation 

ifW(r",r';r) - ^ C„i^„(r", r'; r). (114) 
nez 

where 

C„ = e*2^""' (115) 

and 

li„{t",r';T) 

J-oo \ihsmujTj 

In the coefficients C„ of (jllSp a' can chosen to be an arbitrary real number. In particular, choosing a' = —a 
one observes that the magnetic flux only appears as a pure phase factor exp{— ia(i^" — + 27m)} in (|116p . 
The choice a' = pk — a furthermore shows that the screw dislocation has a similar effect, that is, it only 
appears as a phase exp{— 1(0: — k(5){d" — t?' + 2?™)} in (|116l) . Evidently both belong to ?7(1), the one- 
dimensional unitary representation of ■ni{M.). The result is consistent with the Laidlaw-deWitt-Schulman 
theorem albeit the unitary factor is not unique. 



2T:ih sinwr 



exp 



iMoj 

2h 



(r'^ + r"^) cot LOT 



5 Energy spectrum and wave functions 



All the information concerning the energy spectrum and associated radial wave functions for the bound 
states in two dimensions is contained in the radial propagator (I107p . In order to extract such information 
out of ()107p . we make use of the Hille-Hardy formula (GR: 8.976.1 in [26]). 



2^xyz\ [xyz 



1- z 



where lII^^ (x) is the Laguerre polynomial related to the confluent hypergeometric function as 



exp 



1 , , 1 + z 
'2^" + ^)T3^ 





x + y' 


— exp 


2 



E 



7 1 r" 



^ r(n + /i+ 1 

n— ^ 



Ll^\x)Li^Hy) (117) 



r(n + A* + l) 



F{—n, fi + l; x). 



Notice that the confluent hypergeometric function is in general an infinite series, 

, , Tib) ^ T{a + s) , 
F a, b; x) = > I \ x' 

1 (a) 1 (0 + n) s! 



(118) 



(119) 
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which converges only for \x\ < 1 and becomes a polynomial for any x when a — 0,-1,-2,.... Letting 
X = {Muj/ny^, y = (Afcj/n)r" 2, and z = g-^*"^^ in (fTTT]) . we write (fTOT]) as 



i?,„(r",r';T) = -^f—rV'j cxp 



With this radial propagator the /c-propagator ()llQp can be cast into the form 

oo oo 

K('Hr",e";r\0';r)= ^ ^ ^™„(r", ^?")V4J^', e-^^^""/^ (121) 

m— — oo n— 

where 

E,nn = ?iw(2n + pi{m) + 1) + m^iw (122) 

and 

Substitution of (PT|) into (|122p results in the energy spectrum for the particle bound in two dimensions 

Ernn = ^"^u + 1 + ^/ 4:{m + a- /3ky + a'^ - 1 + + mhu} (124) 

where n = 0,1,2,..., m = 0,±1,±2, Adding the continuous spectrum (155]) as well as the ignored Vb 

yields the full spectrum of the system, 

E„ink = + + {Pk - a)huj. (125) 

6 Concluding Remarks 

In concluding the present paper, we would like to make some remarks on the discrete energy spectrum (I124p 
for the two-dimensional motion around the dispiration. First we examine special cases. 

(i) The Landau levels: The presence of the uniform magnetic field is unimportant for the study on the 
dispiration. However, if there are no dislocation, no disclination, no magnetic tube, no external short-ranged 
repulsive and long-ranged attractive potential but the uniform magnetic field, that is, if a = /3 = k = 0, 
a = 1 and ujq = 0, then, as is expected, we have the Landau levels, 

= 2huj^ j^fi + _ j + __ (126) 

where eB/{2M) and n^n + (|m| +m)/2 = 0, 1,2, • • •. 

(ii) The screw dislocation spectrum: In the absence of the uniform magnetic field, the disclination and 
the inverse-square potential, i.e., in the case of w = 0, cr = 1 (i.e. k = 0), the discrete energy spectrum for 
the two-dimensional motion (with fixed k) becomes 

E„in = hujo C^n + 1 + \m + a ~ l3k\) , (127) 

which shows that the effect of the Burgers vector b — 2ttP is practically identical to that of the magnetic 
tube (the Aharonov-Bohm effect) as pointed out in ref. 11 . In comparison with Wilczek's anyon model |33j . 
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the screw dislocation plays a similar role of an anyon by generating a fractional spin. From the geometrical 
point of view, as we have seen earlier, the line dislocation causes torsion. This spectrum explicitly shows 
that a source of torsion generates a spin effect [31] [32] . In the limiting case of the vanishing flux tube and 
the diminishing Burgers vector, we have simply the harmonic oscillator spectrum, 

En = noJo{n + l), (128) 

where n = 2n -\- \m\ = 0, 1, 2, • • •. The harmonic oscillator potential is important in the present dislocation 
spectrum, without which two-dimensional bound states cannot be formed in the vicinity of the defect. In 
the calculation of the partition function for an anyon gas, the oscillator potential has played a role of the 
regulator for taming divergences [33| [S^ • 

(iii) The wedge disclination spectrum: If a = 0, /3 = 0, and Q = 0, then the two-dimensional discrete 
spectrum takes the form, 

-!-V4m2-hcr2-H-Kl, (129) 

2(7 J 

which belongs to a particle bound near the disclination by the harmonic oscillator potential plus a repulsive 
inverse-square potential. From (fTS]) it is clear that a = 1 implies the vanishing deficit angle 7 = 0. This 
corresponds to the absence of disclination. If < cr < 1, then 27r > 7 > 0, that is, the medium carries a 
positive curvature at the center of disclination. In comparison with the assumed square inverse repulsive 
potential term with k inside the square root of the spectrum, we see that the negative term with — 1 
represents the effect of an attractive force; we may argue that the disclination effectively generates a short- 
ranged attractive force around it. If cr > 1, the medium is negatively curved, and a repulsive force is created 
around the center of disclination (the saddle point). 

The Schrodinger equations for the harmonic oscillator interacting separately with dislocation and discli- 
nation have been solved by Furtado and Moraes [36 . In the absence of the flux tube a — 0, our dislocation 
spectrum (|127l) coincides with their result obtain for the dislocation. However, our disclination spectrum 
(|129p with K = is not in agreement with theirs. 

Finally, we examine for the case of disclination the difference between the result from the Schrodinger 
equation and that of path integration. Since some errors are involved in [36j , we present our own solution of 
the Schrodinger equation which is slightly different from that in [35]. The line element ([T]), if /3 = 0, reads 

ds^ =dr^ +(7^r^d0^ +dz^. (130) 

Although the Laplace-Beltrami operator can be used to write down the Schrodinger equation as in 1361 , we 
choose to let cj) = (t9, and express p30p as 

ds^ ^ dr^ + d<j)^ + dz^ , (131) 

which is identical with the fiat space line element except for < 4> < 2T:a. Then it is rather trivial to write 
down the Schrodinger equation in terms of coordinates (r, 0, z); namely 

- 1? (4) ^^w-^i} * + ''<^'* ^ '"'^ 

which can be easily solved for a two-dimensional harmonic oscillator potential V{r) = ^Aloj^r^ with the 
periodic condition ■(/'(r, 2'Ka + (I), z) = ip{r, (j), z). The normalizable solution is obtained in the form, 

V'(r, (l),z)=N ^^^^^^^•i'l-y^-Mu.r^l-^^ ^1^1/'^ F{~nr, 1 + |m|/CT; (A/w/?i)r2), (133) 

with the condition 

1 / W\ E Ph^\ , , 

1 + — -1— + 7771- («reNo), (134) 



2 V a huj 2M 
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which yields the energy spectrum, 



r Iml 1 k'^h 

£^„.mfc = ?iw<^2n, + l + ^^ + ^. (135) 

Here to € Z as follows from the periodic condition. It is apparent that the above spectrum differs from the 
path integral result (jl29p . The term cr^ — 1 inside the square root of (|129p lacks in ()135p . Furthermore, 
the wave functions are different. The wave functions (jl33|) with m 7^ vanish at r = 0, but the function 
with m — Q remains to be non-zero. This is in contrast with the fact that the radial wave function (|123l) 
obtained by path integration vanishes at r = for all values of m. In this treatment the nonvanishing 
singular curvature at the disclination center r = plays no role. The Schrodinger equation may have to be 
modified so as to accommodate the curvature effect. 

In general, the Schrodinger equation in curved space is written in the form, 

\2M^ + T^c(r) + F(r) j V(r, t) = i?i^^(r, t) (136) 

where A is the Laplace-Beltrami operator, Vc{r) is the potential due to the curvature effect and V^(r) is any 
external potential. Historically, Podolsky [37] defined the Schrodinger equation in curved space without the 
curvature term, i.e., Vc{v) = 0, but, as Schulman fl4 puts it, there is no reason, other than the prejudice 
for simplicity, to ignore the curvature term. Comparing with Feynman's path integral, DeWitt |38j has 
proposed that the scalar curvature term is needed, that is, Vc(r) — gfi'Riv) where 5 is a constant. More 
recently, however, Kleinert 39| has argued by using a quantum equivalence principle that there is no need 
of the curvature term. On the other hand, viewing that the motion in a two-dimensional curved space as 
a constrained motion on the curved surface imbedded in a three-dimensional Euclidean space, Jensen and 
Koppe [10] , da Costa [JT] and others have argued that the Schrodinger equation on a curved surface carries 
in it an effect potential due to the Gaussian curvature K(v) and the mean curvature H{r). 

Our path integral calculation necessitates an effective potential of the inverse square form which is due 
to neither the scalar curvature nor the Gaussian curvature. In a forth coming paper |42| it will be shown 
that the path integration in a conical space with K — and H = y/1 ~ / (2(7r) for r 7^ is compatible 
with the Schrodinger equation modified with the mean curvature H of the conical surface as suggested in 

m- 
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Figures 




Figure 1: Make a thin cylindrical tube along the central line and remove a wedge of angle 7 > 
from the main body. 




Figure 2: Close the open hps in such a way that a screw dislocation is created along the center 
line. The resulting medium is what we call the field of a dispiration. 
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